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ABSTRACT 


A  family  of  generalized  models  are  presented  for  analyzing  data  or 
systems  with  stochastic  properties.  Previous  work  in  this  field  is  pre¬ 
sented  as  well  as  a  new  finite  range  frequency  distribution  function  with 
nomograms,  tables,  examples,  and  recent  extensions  to  a  generalized  family 
of  methods  and  models.  The  analytical  techniques  are  related  to  statistical 
mechanics  and  were  developed  specifically  for  analyses  of  weapon  systems. 
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SECTION  I 


INTRODUCTION 

In  solving  problems  involving  stochastic  processes,  the  systems  analyst 
and  engineer  often  have  to  choose  a  probability  distribution  function  that 
will  fit  the  experimental  data.  One  method  is  to  use  one  of  the  typical 
"text-book"  distributions.  These  are  generally  tabulated,  with  established 
computational  procedures  for  finding  moments.  This  approach  has  proven  useful, 
especially  when  higher  moments  of  the  distribution  are  needed.  However  in 
many  cases,  (a)  only  the  first  moment,  or  the  first  and  second  moments,  are 
needed  in  the  analysis,  (b)  there  may  be  obvious  similarities  to  given  data 
and  curves  encountered  in  previous  situations,  and  (c)  there  may  be  an  ad¬ 
vantage  in  using  a  standard  distribution  that  will  allow  ease  in  some 
predetermined  calculations  and  tests.  In  such  cases,  it  may  be  desirable  to 
devise  distribution  functions  that  fit  the  analysis,  rather  than  encumbering 
the  analysis  with  unnecessarily  complicated  functions  and  fitting  procedures. 

In  this  report,  a  class  of  distribution  functions  is  derived  and  it 
is  shown  that  the  frequency  function  developed  by  Braswell  and  extended  by 
Manders,  References  1  and  2  (see  the  Appendix),  is  a  particular  member  of  a 
family  of  these  distribution  functions.  Procedures  and  methods  of  finding 
parameters  of  these  distributions  are  presented.  To  show  the  simplicity  and 
usefulness  of  this  method,  the  parameters  of  the  distribution  of  Braswell  and 
Manders  are  re-derived,  using  the  more  general  techniques.  Finally  an 
example  of  application  of  these  types  of  distributions  is  given. 
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SECTION  II 


CLASS  OF  FREQUENCY  FUNCTIONS 


Suppose  there  is  given  a  function  h(x),  such  that 
dh(x) 


dx 


>  0 


for  Xq  £  x  <_  x^ 


Then  the  doubiv  truncated  cumulative  distribution  function 


F(x)  =  < 


f  0  for  x  <  Xjj 

h (x)  -  h(xQ) 
h (x^)  -  h(xQ) 

1  for  x^  <_  x 


for  x.  <  x  <  x, 

0  1 


is  a  proper  cummulative  distribution  function  (cdf)  with  finite  domain 
associated  with  positive  probability. 

The  corresponding  probability  distribution  function  is  given  by 


f(x) 


0  for  x  i  [x0,x1] 

L  h(x1)'-3h(x0j  for  X  £  [x0,xl] 


The  mean,  p ,  is  given  by 


u  = 


xf(x)  dx 


=  xF(x) 


X1  fXl 


F(x)  dx 


x0  xo 


Since  F(Xq)  =  0  and  F(x^)  =  1, 


,-r 


U  “  x,  -  I  F(x)  dx. 

x0 


2 


n 


Similarly,  the  second  moment,  u2»  is  given  by, 


W,  =  x  F(x)  -  2  xF(x)  dx 


f1 

>  x. 


Since 


d  fx 

^  I  F(t)  dt  =  F(x),  Equation  7  may  be  written, 

J  V 

0 

x 


P2  =  x  F(x) 


1  fx 

-2x 

e  v 


F(t)  dt 


X,  X.  X 

xo  xo  xo 


F(t)  dtd< 


Therefore,  a  useful  form  for  y ^  if  F(x)  is  twice  integrable  is 


U2  =  xx  -2xx 


1  f  1  fX 

F(t) dt  +2  I  F(t)  dtdx 

*0  xo  xo 


Using  Equations  6  and  9,  y2  nia y  be  expressed  by 


X j  x 

-  x  y  -x.J  F(t)dt  +  2  [  F  (t)dtdx 


x0  X0 
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The  absolute  deviation,  y,  is  given  by 

x. 


[1 

|x  -  u! 

^  V 


f(x)  dx 


11 


For  Xq  >_  0, 


■  V 

v 

r  1 

xf(x)  dx  +  y 

L 

f(x)  dx  + 

V 

xf  (x)  dx  -  y 

1 1  ^ 

f  (x)  dx  12 


=  -  xF(x) 


+  |  F(x)  dx  +  u[F(u)  -  F(xq)]  +  xF(x) 

X0  xo 


F(x)  dx  -  p[F(x1)-F(y)] 
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r 

j  V. 


F(x)  dx  - 


f  X1 


F(x)  dx  -  u  x  x, 


14 


3 


To  show  che  simplicity  of  use,  the  form  and  parameters  of  the 
frequency  function  of  Braswell  and  Manders  are  derived. 

deax 

Here,  h(x)  =  -  for  x.  =  0  and  x,  =  1 

1  i  U  aX  0  1 

1  +  be 

and,  thus  h'  (x)  >  for  a  >  0,  x  e  [0,  1]  . 

Also,  h(xQ)  =  =  x~+  b 

h  (x.  )  =  h(l)  =  — — — - 
1  1  +  bea 

Therefore , 

H(x)  for  0  <_  x  <  1 

To  get  F (x)  into  the  form  previously  used,  we  examine  H(x). 


H(x)  =  IefL-_lI(J-L±e!) 

(J  +  beax)(ea  -  1) 


[1  +  e~*(1+2s)]  fj  -  e  2UX] 

[  1  -  e_2a  )  [1  +  ea(1'2x-26>'] 


for  0  <  x  <  1 


,  a  ,  t  a  +  2  in  b 

where  a  =  -  v  and  5  =  - - - 

l  za 


a(l-26) 

Multiplying  Equation  17  by  VIY 

Ot  ^  1  —  L  o  ) 
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fl  ^e-a(l^)Hea(l-26)  _  ea(l-2x-26)] 

W  '  ,  a(l-25)  -a (1+25)  ,  .  a(l-2x-26)  . 

I  e  -  e  ] [ 1+e  ] 


for  0  <  x  <  1 


which  is  the  result  obtained  by  Braswell  and  Manders. 
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The  mean,  u,  is  given  by  Equation  6,  therefore 


1  rl 


u  -  xF(x) 


F(x)  dx 


Thus,  y  may  be  written  as 


V  *  1 


/I  ±  bea\  f1  1  -  eaX 
\1  -  ea  /•'01  + 


The  integral  in  Equation  20  is  evaluated  in  Reference  1. 

Due  to  desired  generalizations  for  the  cdf  of  Braswell  and  Manders 
presented  later  in  this  paper,  it  is  desirable  to  be  able  to  evaluate  integrals 


of  the  form 


/  1  cx , 
e  dx 

I  ..  aX 

x0  b+e 


Thus,  the  transformation,  n  “  e  gives 
ax-, 

.  re  c/a-1 , 

i  *  If  n - 4a 

a  I  ex_  b+n 


axn 

For  b  >  -  e  u  , 


_ 1  1  1  1  n 

n1_c/a(b  +  n)  ’  b  n1_c/a  "  b  {b  +  n) 


Therefore , 


ax  j 

e  c/a. 


,  e  c/a, 

dx  _  1_  x  dx 

1-c/ a  ab  J  ax  (b  +  x) 

IF  '  /-» 


Evaluating  the  first  integral. 


I  *  —  (  -Jl 
1  ab  j  ax,,  1-c 


1  c/a 
x 

be 


ax0  for  c  /  0 


ax0  for  c  =  0 


5 


w 


There fo re ; 


X1 


i  a*i  axo 

££  (e  -  e  ]  for  c  ^  0 


b  [xi  _  xo] 


for  c  =  0 


26 


The  second  integral  to  be  considered  is, 
ax^ 

-  i_  f  3  xc/a 

l2  ab  J  aXo  (b+x)  dx 

e 


27 


If  c  =  0,  then 


i  =  —  |  4*. 

2  ab  I  axn  b+x 


28 


*2  =  ^b  '’n 


ax^ 


axr 


29 


If  c/a  =  n,  some  positive  integer,  then  I-^is  given  by  Equation  25,  and 


axi 


2  ab 


dx 


30 


ax^ 


-  5b  J 

r  e  +b 

(x-b)n 

ax0 

e  +b 

X 

1 

axi 

r  e  +b 

y11  (nl 

ab  ' 

axo 

e  +b 

dx 


31 


!  =  £*> 
2  ab 


axi 

n  r  e  +b 


K0fb 


e  +b 


x^1 


ax0 
e  +b 


dx  32 


6 


7 


8 


the  fewer  number  of  terms  necessary  in  the  evaluation). 

An  immediate  result  of  the  foregoing  discussion  is  a  method  of 
evaluating  slightly  more  general  integrals,  those  of  the  form 


Here 


I 


1 

a 


e  c/a 

x _ 

ex0  x(b+x)a 


dx 


37 


38 


1 

a 


X1 

r 

)  xo 


+b 

+b 


(x-b)^3-1 

a 

x 


dx 


39 


yc/a-l\  eXl+b 

-  i  n  n(  )  (-b)^3"1^  dx  40 

a  \  ,  /  J  x0 

^  e  +b 

Again,  if  c/a  is  a  positive  integer,  there  will  be  only  c/a  +  1  terms 
in  Equation  40. 

It  is  of  interest  to  note  that  the  foregoing  integrals  are  generaliza¬ 
tions  of  many  of  the  tabulated  integrals.  Therefore , the  evaluation  of 
these  integrals  is  of  some  importance  in  itself. 

Returning  to  the  determination  of  the  mean  and  absolute  deviation  of 
rbe  Hisfri’Kii.for,  oc  ”raswe11  and  Manders,  the  mean  is  obtained  using  Equa¬ 
tions  20,  2b,  and  31.  The  absolute  deviation  is  obtained  by  use  of  Equa¬ 
tions  14,  2b,  and  31.  When  a  =  -  -■  and  5  =  2a^°^  *  t*le  results  in 

Reference  l  are  obtained. 

For  completeness,  the  frequency  function  of  Braswell  and  Manders  should 
be  examined  for  the  limiting  cases: 

(a)  Lim  F(x)  (b)  him  F(x) 

a-0+  b-*~l+ 


8 


Q 


For  case  (a) 


Lim  F(x)  *  Lim 


a+0, 


a+0, 


/.“-ll 

/l+bea  1 

.ea-l  i 

il+beaxl 

If  b  4  -1,  then 


Lim  F(x)  «  Lim 


a+0. 


a+0. 


eax-l 

ea-l 


Therefore,  by  L'Hospital's  rule 

Lim  F(x)  =■  x 
a+0+ 

Thus  F(x)  tends  to  the  cdf  of  a  random  variable  uniformly 
distributed  in  the  interval  {0,1}. 

In  case  (b) 


41 


42 


43 


or 


Lim  F (x) 
b+-l 


Lim  l+bea 

b+-l,  ...  ax 
+\l+be 


Lim 
b-*— 1 . 


F(x) 


1 


Thus  x  is  a  variable  associated  with  a  deterministic  process. 


44 


45 


A  useful  extension  of  the  frequency  function  of  Eraswell  and  Manders 
given  in  Reference  1  is  given  by 


h(x)  - 


cx 


b  +  e 


ax 


46 


with  x  =  0,  x,  =  1,  b  >  -l. 
o  ’1  ’ 


For  monotonicity  of  F(x),  must  be  greater  than  zero, 


V"  x  for  x  e  [0,1] . 


Since 


dh(x)  m  £e 


cx 


dx 


ax 


cx  ax 
ae  e 

(b+eaX) 2 


the  condition 


b+e 

>  0  will  hold  if 
dx 


b+e 


ax 


x  for  x  r  [0,1] 


47 


9 


10 


ax 

£ 

Since  w(x)  =  — — —  is  monotone  increasing  for  x  £  [0,1],  the  maximum 
b+eax 

of  w(x)  occurs  at  x  =  1,  and  therefore,  monotonicity  is  insured  if 
a 

c  >  - — —  ,  (b  >  0  and  c  >  a  will  insure  monotonicitv). 
b+e3 

It  should  be  noted  that  h(x)  given  in  Fquation  46  is  a  generalization 
of  the  energv  distribution  functions  of  statistical  mechanics  as  given  in 
Reference  3.  The  Maxvell-Roltztnann  distribution  is  of  the  form. 


f(x) 


x  x/kt 
e  e 


The  Bose-Einstein  distribution  is  of  the  form 


f(x)  = 


a  x/kt  . 
e  e  -1 


48 


49 


and  finally,  the  Fermi-Dirac  distribution  is  given  by 


f(x) 


a  x/kt. 
e  e  +1 


50 


Using  h(x)  given  in  Equation  49,  the  cdf  takes  the  form 
I’  0  for  x  <  0 


F(x) 


cx  l 

e 

1 

1 

1 

k  ,  ,  ax 
l  b+e  j 

b+1 

I  c  ^ 

_L_ 

\b+e  ;  l b+1 

for  0  <  x  <  1 


51 


'  1  for  x  >  1 

The  function  given  in  Equation  15  is  a  particular  example  of  a  family 


of  functions  given  by 
cx 


h(x) 


de 

,,,  ax.  a 
(b+e  ) 


for  x^  =  0  and  x^  »  1 


52 


Using  this  function,  the  corresponding  cdf  is  given  by 


10 


To  evaluate  the  parameters  of  this  distribution, 


integrals  of  the  form 


54 


must  be  evaluated  and  the  evaluation  of  these  integrals  is  greatly  simplified 
(thus  increasing  the  usefulness  of  the  distribution)  if  c/a  =  n,  a  positive 
integer  and  a  *  k,  some  positive  integer. 


In  this  case.  Equation  54  becomes 


55 

56 

57 


There  are  four  parameters  associated  with  this  distribution:  a, 
b,c,  and  k  (n  is  determined  by  a  and  c)  and  formulae  have  been  developed 
earlier  for  evaluating  the  mean  and  absolute  deviation.  Therefore,  if 
any  two  of  the  parameters  are  fixed,  the  mean  and  absolute  deviation 
can  be  used  to  determine  the  other  two.  By  doing  this,  six  families  of 
distribution  functions  are  generated. 


11 


For  example 

(  0 


if  b 


1 ,  then 


=  1,  a  -  1,  xQ  -  0,  xx  - 


F(x) 


for  x  <  0 


I  enx 

1 

1 (l+ex)k 

2k 

n 

\ 

e 

1 

i(l+e)k 

2k] 

^  1  for  x  >  0 


for  0  <  x  <  1 


58 
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SECTION  III 


RANDOM  NUMBERS 

Often,  it  is  necessary  to  use  random  numbers  selected  from  a  distribution. 
The  cdf  of  Braswell  and  Manders  may  be  easily  inverted  to  give  random  numbers  from 
that  distribution  if  random  variables  uniformly  distributed  in  {0,1} 
are  available. 

If  F(x)  is  monotone  increasing  cdf,  and  r  is  a  variable  from  a 
uniform  {0,1}  distribution,  then 
y  =  F_1(x) 

is  a  random  number  from  F(x). 

For  the  FRPDF  of  Braswell  and  Manders,  the  cdf  may  be  inverted  in  closed 
form,  thus  making  generation  of  corresponding  random  numbers  easy. 
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TABLE  I. 


RANDOM  NUMBERS  USING  THE  FRPDF  OF  BRASWELL  AND  MANDERS 


0.  24158 

0.60196 

0.25900 

0,97604  ! 

0.81010 

0.32882 

0.54244 

0.02165 

0.99494  i 

0.35506 

0.57602 

0.21315 

0.83320 

0.37908 

0.46429 

0.46475 

0.47708 

0.07142 

0.19321 

0.46652 

0.59312 

0.59904 

0.55056 

0.37447 

0.75369 

0.74804 

0.00798 

0.53492 

0.38361 

0.84148 

0. 19745 

0. °6575 

0.11968 

0.55770 

0.64208 

0.42256 

0.74125 

0.65361 

0.77012 

0.59652 

0.03307 

0.38003 

0.21965 

0.22837 

0.02158 

0.32717 

0.51644 

0.43901 

0.31461 

0.38016 

0.73331 

0.38885 

0.23166 

0.54247 

0.09798 

0.58225 

0.18766 

0.31569 

0.79382 

0.04368 

0.26919 

0.25216 

0.43510 

0.26635 

0.00619 

0.50120 

0.67743 

0.35515 

0.32886 

0.48071 

0.96026 

0.04537 

0.41694 

0.79842 

0.68748 

0.37894 

0.50911 

0.43165 

0.87127 

0.67873 

0.99545 

0.08539 

0.21373 

0.58254 

0.36961 

0.74520 

0.19261 

0.47631 

0.79148 

0.19704 

0.19635 

0.29674 

0.40224 

0.36312 

0.51801 

0.75558 

0.82781 

0.99734 

0.18026 

0.00213 

CA=0. 10000,  B^O. 90000) 


Example ; 

It  is  illustrative  to  show  how  the  generalized  class  of  frequency  functions 
can  be  used  in  systems  analysis. 

A  sample  of  times  to  complete  jobs  on  a  computer  was  taken  (Table  II). 

A  useable  empirical  distribution  was  needed  to  fit  such  samples. 

To  handle  this  problem,  the  probability  distribution  associated  with 
the  generalized  cdf  vas  analyzed. 

It  was  noted  that  these  jobs  fell  into  two  categories.  The  first 
category  was  short  jobs,  i.e.,  jobs  that  took  less  than  one  minute  to 
complete.  The  other  category  was  long  jobs,  i.e.,  those  that  used  more 
than  one  minute  to  complete. 
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The  shorter  jobs  far  outnumbered  the  longer  jobs  at  the  installation 


from  which  the  data  was  taken.  Also,  it  was  known  that  the  longer  jobs 
were  usually  "production  runs"  of  tested  programs,  so  that  they  could 
be  handled  differently  than  the  short  test  jobs. 

Thus  the  data  of  interest  concerned  the  shorter  running  programs. 
The  associated  pdf,  f(x)  may  be  written, 


f  <x)  = 


,,  .  ax.k 
(b+e  ) 


kae 


b+e 


(b+ea)k 


(b+1) 


59 


Since 


ax 

6  k.36 

is  an  increasing  function  of  x,  and  c- 


...  ax.k 
(b+e  ) 


b+e 


ax 


is  a 


decreasing  function  of  x,  it  was  seen  that  the  important  middle  range 
values  of  the  sample  could  be  fitted  by  varying  some  of  the  parameters. 

Since,  only  means  and  absolute  deviations  were  easily  computed, 
the  cdf  given  by  Equation  58  was  adopted.  Then 


r 

1 

nx 

e 

keX 

|^(l+eX)k_j 

n  — 

l+eX 

n 

e 

i 

_(l+e) 

k  ' 

2k  . 

A  table  of  means  and  absolute  deviations  for  various  n  and  k  is  shown 
in  Table  II  and  used  to  pick  a  value  of  n  and  k  to  agree  with  the  estimated 
mean  and  absolute  deviation  of  the  sample,  which  were  u  =  .2534,  and  y  =  .17. 
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TABLE  II. 


1 


U  AND  y  AS  A  FUNCTION  OF  n  AND  k 


K 

n 

=  1 

k 

n  = 

2 

u 

_Y _ U 

_ 

L 

u 

y 

0 

0.577 

0.235 

0.652 

0.219 

1 

!  0.475 

0.242 

0.235 

2 

0.640 

0.201 

0.506 

0.240 

3 

0.494 

0.238 

4 

0.429 

0.227 

0.624 

0.200 

5 

0.378 

0.217 

0.497 

0.229 

6 

0.333 

0.206 

0.428 

0.225 

7 

0.296 

0.188 

0.374 

0.217 

8 

0.263 

0.174 

8 

0.329 

0.196 

9 

0.236 

0.156 

0.291 

0.189 

n 

=  3 

n  = 

4 

0 

0.174 

mm 

0.764 

m 

1 

0.669 

0.729 

2 

0.613 

0.688 

1  3 

0.536 

0.635 

!  4 

0.383 

0.566 

5 

5 

0.457 

0.228 

1  6  ! 

0.608 

6 

7 

0.494 

7 

8 

0.423 

8 

0.592 

0.204 

9 

0.369 

9 

0.487 

0.220 

n 

=  5 

n  3 

6 

i  0 

0.802 

■ 

0.831 

mssmm 

1 

0.776 

fXSSSi 

0.811 

i  2  ; 

0.744 

mSEm. 

0.787 

Hi.  |  ■ 

3  i 

0.705 

0.758 

■  si  •  •  • . 

4  1 

0.657 

0.722 

5  ; 

0.594 

BrAM 

I 

5 

0.678 

9>'i  -  ■ 

6  1 

0.504 

0.233 

6 

0.620 

■SSl  - 

■ 

7 

0.548 

1 

I 

8 

0.420 

0.207 

JL. 

■ 
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For  the  given  sample,  n  =  1  and  k  =  8  were  chosen.  To  test  the 

goodness- of- fit  of  the  fitted  distribution,  the  data  were  divided 

2 

in  five  group  and  an  X  test  performed.  See  raw  data  in  Tables  III  and  IV. 

The  regions  considered  were  (data  in  Table  V) . 

I:  0  <_  x  <  0.1 

II :  0. 1  <_  x  <  0.2 

III:  0.2  <_  x  <•  0.3 

IV:  0.3  <_  x  <  0.4 

V:  0.4  <  x  <  1.0 
2  (y  -g 

Setting  X  _  L  Jm  6m  ,  where  y_  =  number  of  data  points  in  interval 
m=l  „  m 

sm 

m,  and  g^  =  expected  number  of  points  in  interval  m.  g^  was  found  using 

the  tabulated  value  of  F(x) ,  and  x  was  found  to  be  4.087.  Since  there 

were  five  intervals,  and  two  estimated  parameters,  the  variable  x2  should 

belong  to  a  chi-squared  distribution  with  two  degrees  of  freedom.  The 

value  of  4.087  indicates  a  significance  level  of  about  85  percent. 

Much  higher  values  of  x  are  found  if  a  normal  distribution  is 

hypothesized,  due  to  the  fact  that  o2  *■  .044.  This  small  value  of  o? 

causes  g^  and  to  be  small,  thus  increasing  the  value  of  x  •  The  given 
2 

X  for  this  problem  could  be  greatly  reduced  if  the  restriction  of  integral 
n  and  k  were  removed,  or  if  larger  values  of  integral  n  and  k  were  considered. 
CONCLUSION 

The  scheme  presented  in  this  paper  may  be  extended  to  use  for  different 
"basis"  functions  for  distributions.  To  use  this  method  the  analyst  should: 
(a)  Determine  if  there  is  something  to  be  gained  by  using  this 

technique.  If  the  problem  under  consideration  is  not  made  easier, 
he  should  not  pursue  this  technique. 
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TABLE  III.  SEQUENCE  OF  COLLECTED  DATA 


1 

Observation 

Time 

. 

Observation 

_ _ 

Time 

— 

i  Observation 

Time 

Observation 

Time 

1 

0.46 

12 

0.05 

23 

m 

0.04 

2 

0.23 

13 

0.22 

24 

mm 

0.08 

3 

0.35 

14 

0.39 

25 

KH 

4 

15 

0.13 

26 

Hi 

37 

0.53  | 

5 

0.37 

16 

0.88 

27 

0.02 

38 

0.04  i 

6 

0.10 

17 

0.31 

28 

0.08 

39 

0.40  I 

7 

0.07 

18 

0.08 

29 

0.10 

40 

0.13  j 

8 

0.34 

19 

1.00 

30 

0.05 

41 

9 

0.09 

0.49 

31 

0.47 

42 

0.26  j 

IHH 

0.59 

21 

0.15 

32 

0.39 

43 

0.  30  j 

ii 

0.32 

22 

0.27 

33 

i 

0.12 

44 

0.02 

TABLE  IV.  ORDERED  DATA 


TABLE  V.  CLUSTERED  DATA 


Interval 

No.  of  Observations 

0  to  0.1 

13 

0.1  to  0.2 

8 

0.2  to  0.3 

6 

0.3  to  0.4 

9 

0.4  to  1.0 

8 

Parameter  estimates:  u  =  0.25,  y  -  0.17,  a  »  0.044 
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(b)  Pick  a  family  of  functions  that  may  be  parameterized,  and  that 
he  can  manipulate  in  some  way. 

(c)  Determine  the  forms  of  the  associated  cdf's  and  pdf's.  Of  parti¬ 
cular  interest  would  be  critical  points  of  these  functions. 

(d)  Establish  a  method  (tables  of  means,  deviations,  etc.;  explicit 
formulas;  etc.)  of  estimating  the  parameters  of  these  distributions. 

Once  he  has  these  tools  he  may  then: 

(e)  Proceed  to  solve  the  problem  at  hand. 

Objections  to  this  type  of  use  may  be  raised,  and  the  analyst  would 
have  to  determine  the  appropriateness  of  applying  this  method  in  a  particular 
application.  One  point  to  be  considered  in  this  determination,  is  that,  in 
general  problems  that  are  presented  for  solution  are  so  large  or  complicated 
that  a  digital  computer  must  be  used.  If  a  numerical  approach  is  used  for 
parts  of  a  problem,  then  the  analyst  should  be  able  to  use  this  method  if  it 
makes  some  parts  of  the  problem  more  manageable.  But,  since  this  is  the 
only  case  in  which  this  method  should  be  considered,  most  objections  to  the 
technique  can  be  overcome. 
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(The  reverse  of  this  pace  is  Man*' 


APPENDIX 


Much  of  the  information  in  this  report  is  based  on  mathematical 
equations  presented  in  References  1  and  2.  Those  papers  were  published 
only  in  Japan  and  are  not  readily  available  in  the  United  States.  There¬ 
fore,  by  special  arrangement  with  the  publisher,  the  papers  are  included 
here  as  an  appendix. 
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A  NEW  FINITE  RANGE  PROBABILITY  DISTRIBUTION  FUNCTION 
(FRPDF)  WITH  PARAMETERS-NOMOGRAM  AND  TABLES* 


By  Robert  N.  Braswii t 

University  of  Florida 
Gainesville,  Florida 

Clara  Fu-Mei  Mandlrs 
Radiation,  Incorporated  Melbourne,  Florida 


The  objective  of  this  paper  is  to  introduce  a  flexible  finite  range  probability 
distribution  function,  FRPDF.  The  emphasis  is  on  flexibility  of  application 
relative  to  its  simplicity  of  use  and  its  ability  to  fit  varied  experimental  data 
clusters. 

The  FRPDF,  /(*),  (F^x^t,  is  strictly  unimodal  as  shown  in  Figures  2 
through  7.  When  6  =  0,  the  probability  density  function  is  always  symmetrical 
with  respect  to  the  vertical  line  x  =  0.5,  and  the  distribution  function  is  a 
family  of  S-shaped  curves.  When  S> 0,  the  peak  of  the  probability  density 
function  shifts  toward  the  x=l  line  and  when  <5> 0,  it  shifts  to  the  x-0  line. 
When  a  =  0,  the  probability  density  function  is  identically  uniform.  For  the 
large  value  of  «,  a  =  10  or  more  the  peak  of  the  probability  density  function 
becomes  narrower  and  higher  at  appropriate  values  of  x  for  given  values  of 
S.  Then,  6  can  be  considered  as  the  location  parameter,  while  a  can  be 
considered  as  the  shaping  parameter  and  these  parameters  will  give  the  desir¬ 
able  feature;  namely,  flexibility  to  the  probability  density  function. 


1.  Introduction 

The  purpose  of  developing  this  new  distribution  function  is  to  fulfill  the  need  for  a 
flexible  and  easy-to-use  finite  range  probability  distribution  function  for  many  experimental 
problems.  Knowledge  of  the  probability  distribution  of  a  random  variable  is  required  be¬ 
fore  statistical  inferences  can  be  made.  All  real  random  variables,  by  their  very  nature, 
have  a  finite  range.  The  function  used  as  a  starting  point  in  this  paper  was  originally  used 
by  Pearl  to  approximate  the  population  growth  characteristic  of  the  United  States  1  .  Then 
BraswelUC  reformulated  it  into  Task  Operating  Characteristics  (TOC)  curves.  This  paper 
transforms  the  TOC  curve  into  hyperbolic  expressions  as  shown  in  Fquation  (4)  and  into 
the  finite  probability  density  function  shown  in  Equation  (5). 

The  developed  FRPDF  is  like  the  beta  distribution  function  in  that  it  contains  two 
parameters,  a  and  d,  and  that  the  admissible  values  of  the  variate  .V  lie  between  zero  and 
one.  Unlike  the  beta  distribution,  it  is  flexible  and  simple  to  apply  as  illustrated  in  this 
paper.  Fn  general,  this  distribution  is  useful  in  any  application  requiring  statistical  analysis 
of  collected  experimental  data.  The  computer  was  used  to  develop  tables,  curves,  and  a 

*  Received  24,  Feb.,  1970. 

This  work  was  supported  in  Part  by  THEMIS  Contract  DA  FI  ('04  MC  0002. 

[Re.pfu.ntzd  ^nam  Rc.p,  Stat.  A ppl.  Rzd>.,  JUSE,  Vet  77,  Sic.  2,  1970, 
by  peAjmiibfon  the  publtbheA.) 
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nomograph  10  relate  the  mean  and  tne  absolute  delation  of  the  sample  with  the  density 
function  parameters,  a  and  3. 


2.  Mathematical  Basis 


Pearl's  function  for  the  population  growth  of  the  United  States  is  defined  as: 

“  I  +«**•  v  0  (‘'me)  (1) 

where  >-=  number  of  4cople;  a,  b,  and  c  are  parameters;  and  e  is  the  base  for  Naperian 
logarithms.  Braswell's  first  paper  changed  Pearl’s  function  into  the  following  form: 

L'+e 

K  )  \_e“  1  -•  -e J  1  ; 

for  0;£r<l,  a  0,  and  -0.5  -<5i_+0.5. 

Since  Pit)  is  not  defined  when  «=0,  it  is  necessary  to  define  this  function  at  this  point 
by  use  of  L  ‘Hospitals’  rule. 

‘f  A ’(a,  3,  r) 

limP(/)=lim  ‘f  =1 

°-“aD(a,S,') 

where 

N(a,  3,  0  =  [l+e  o,l'-'; 

and 

D(a,  3,  = 

For  more  'completeness.  Equation  (2)  becomes 

=  c  1  +  e  a' 1  ' "'He*" - 2i'  - <!“"-■ 2<1] 

Pit)  ior  «>0  (3) 

=  r  for  rr  =  0  0<f<l, 


A'for,  3,  t ) 


for  a>0 


3.  Formulation  of  the  FRPDF 

For  convenience  change  the  notation  of  the  independent  variable  from  t  to  x  and  P  to 
F,  and  let  F(x)  be  a  probability  distribution  function.  When  one  rewrites  these  equations 
more  explicitly,  one  obtains  the  following  result.  For  a>0,  Pit)  becomes: 

‘X‘  1  —g  a,  I  •  WijQl 

and  for  a  =  0  F(x)=x,  for  0<ar<1.0. 

Equation  (3)  can  be  rewritten  as  follows: 

■  =  cosh  [a(d  +  0.5)]  sinhaJC  -  n 

F(x)  sinh «  cosh  [><.c  +  S-0.5X]  a>  (4) 

.  —  x  for  a  =  0  0£a:<l. 

Since  Fix)  has  the  following  properties,  it  is  a  valid  probability  distribution  function-3-: 

(a)  F(0)  =  0,  F(l)=  1 

(b)  It  is  a  nonaccreasing  function  of  x: 

Fix, )L  Fix.  )  for  x  <x2 

(c)  It  is  continuous  from  the  right: 

Fix  )^  Fix) 

i.e.,  if  .ti <.r2,  then  P{x,<  A’  J.v-}  Fix •)- Fir,). 

For  a  continuous  random  variable,  the  probability  density  function  is  F’{x).  For  a>0, 


for  a> 0 
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/(.*)  = 


_  /  o-cosh  n«(3  +  0.5)_]  cosh  aiu- 0.5)  |  \  / 
\  sinh  a  I  \ 

Consequently  the  complete  form  for  the  FRPDF  is: 

=  C("’  ^  cosh2[or(x  +  3-0.5)J  1 

J{X)  =  1 .0  for  a  =  0  OiT.v: '  1 .0 
“0  elsewhere 


=  G(«,  3 ) 


cosh-  |.<i(.v  i  3  -  0.5)  J 


for  «>0  0  I.v:'l 


G(a  8)—  a  cos^  r<«(5  +  0.5)J  cosh  fo(3  -  0.5)j 
sinh  a 

The  useful  parameter  approximate  .is  of  the  FRPDF  are  the  mean  and  absolute  deviation 
of  the  samples.  The  mean  and  the  absolute  deviation  as  a  function  of  a  and  3  is  derived. 
It  will  also  be  shown  that  the  variance  is  finite  and  that  the  moment  generating  function  exists. 


4.  The  Mean 

The  mean,  /i,  is  by  definition: 

/i  =  £(x)  =  |  xf(x)dx. 

Substitute  the  Equation  (5)  into  Equation  (6),  for  a>0,  0<a<1, 

/l(a'  »- !.**«•  cosh2Ca(x+ 3— 0.5)3  d* 

~G{a,  cosh-’ [a(x+  d- 0.5)1  d*' 

Let  ;  =  rt(x  + 3-0.5). 

Then  /t  becomes 

rf«i‘fl. 5)  d~  roito. si  1  j,  T 

/*r,  J)-G(«r,  3)  \  -  +  \  (0.5-3)  ‘  1 

Ua'j-o.51  arcosrr;  a  Ja-s-o.5>  cosh-;  a  J 

J  z  ,  (0.5-3)  fn5*°- 51  1  .1 

’  L  a2  L  cosh-;  z  a  L  j-o.d  cosh2 ;  ^J‘ 

From  integral  tables  [4J, 

"<“■  »-«*•  »  l  [«  “"h WS+0.W+IO*  ] 


for  ar>0  and 


The  mean  is  therefore 


i  =  \  xdx :  =“  0.S  for  «  =0. 


At..  »-)■  ***  Mog  ]  for  «>. 


[  =  0.5  for  «  —  0 

and  G(<>.-,  3)  is  defined  in  Equation  (5). 


5.  Absolute  Deviation  [5] 

The  absolute  deviation,  denoted  r,  is  defined  as 

r(tt.  o)-  \  | .V  / 1 !  Kx)<lx 
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Table  1 

Mean 

as  a  function  of  a 

and  d 

-  .  5 

-0.4 

-0.2 

-3.1 

0.0 

+  0.1 

* 0.2 

+  0.3 

+  0.4 

+0.5 

0.5696 

0. 5570 

0. 5u 36 

0.  52*35 

0.5143 

0.5000 

0.4851 

0.4705 

0.4504 

0.4430 

0.4304 

.• 

0. 6872 

0.0595 

0 . 6  2  4  7 

0. 5 o  7 o 

C . 5445 

0. 5000 

0.4555 

0.6130 

0.3743 

0.3407 

0.3128 

0. / 7 34 

e.i  !  3 

0  .  u  >  .  4 

0.5693 

0.5060 

0.4367 

6.3655 

0. 3052 

0.2615 

0.2264 

4 

6 . 5  2  7  •* 

0. 7cZi 

0.7357 

0  ,  t>  U  4  4 

t. 5345 

0.5060 

0.4155 

U.  3  359 

0. 26G5 

0.2112 

0. 1726 

0 .  c  u  1  5 

0 . 3  2  0  5 

U  .  7  5  .  6 

O.uSO# 

C  .  5  V  2  6 

U.  5')uU 

u .  4u  74 

0.3193 

0.2412 

0.1795 

0.1385 

b 

0 . 3  -  *  5 

0.8-U 

6.  7/  -.U 

6.6;  7 

0 . 5  5  o  5 

0.5^00 

0.4035 

0.  51C3 

U.  2260 

0.1586 

0.1155 

7 

0 . 0  -  1  c 

o .  d  » :•  / 

- .  7-34 

0.4  .h5 

0.53-4 

6 . 5  o  6  0 

0.4616 

0 . 3055 

0.2166 

0.1443 

0.0990 

8 

o .  .  i :  4 

J.6lb'j 

o .  7  -  • : 

6  .  L-  )  7  J 

0.  V')9  J 

0.5 GUO 

0.400/ 

0.3636 

6.2167 

0.1340 

0 . 1)  -  l>6 

y 

o .  «■ :  3  o 

0. 8  7  ;-o 

0.7)59 

0 . 6  j  3  4 

0.54-7 

0.53-0 

0.4003 

0. 3C1G 

0.2070 

0, 12o4 

O.U770 

1  J 

0 . 0  .  7 

0 . 3  /  ‘.  3 

0.7)54 

0.6'j  jl 

0.5749 

0 . 5uC0 

0.4601 

0.3069 

0.206G 

0.1207 

0.0693 

11 

C .  9  3  #  0 

6.8  .i  56 

0.7)73 

0.69)5 

0.5)09 

0.5003 

0.4601 

6.3005 

0.2030 

0. 11G4 

0.0G30 

:  2 

C « 'j  *  J  2 

G.--7J 

u, 7 1  :o 

0.0tj7 

0 . 6  l  L  0 

0. 5000 

0.4600 

0. 3003 

0.2020 

0.1130 

0.0578 

.  \ 

0 . 4  4  a  7 

6  .  3  S  'j  6 

j .  7  j ;  7 

Q.i  >  ) 

0.1.6.4  a 

O.SvUO 

0.4000 

0.3661 

0.2013 

0,1104 

0.0533 

:  4 

0.4505 

0.0517 

0.7)91 

O.vfjJ 

U.bti 0 

U . 5600 

u . 4000 

0. 3001 

0. 2UC9 

0.1083 

0.049S 

. 

0.9538 

9.3933 

C.75a4 

0.76.40 

C  .  v.6  " 

1.5GC0 

Q . 4060 

0.3000 

0.2006 

0.1067 

0.06G2 

.  ■„ 

-.9567 

0.354O 

«‘,7  1 0 

0. /u-0 

O.uf 

6 .  :o  6  0 

0.6006 

0.30U0 

0.2004 

0.1054 

6.0433 

.  7 

0.9592 

0.8*):»7 

0.7797 

0.  /*•«.) 

o .  u  u  o ; 

U.5'.<3  0 

0.6090 

0.3000 

0.2003 

0.1063 

0.0403 

.9 

f<  •  5  G 1 5 

0.8 JuS 

C. 7558 

0. 7f vO 

0.60C* 

0, 5CC0 

C.4C0O 

0.3000 

0.2002 

0.1035 

0.0335 

19 

0,9635 

0.8972 

0.7399 

O.-’w^ 

9,6000 

0.5000 

0.4C00 

0.3000 

0.2001 

0.1028 

0.0365 

JU 

0.9653 

0.8977 

0.79)9 

Q.  70 

.  6000 

0.5000 

0.4000 

0.3000 

0.2001 

0.1023 

0.0347 

Table  2 

Absolute  deviation 

as  a  function  of  a  and  ) 

-0.5 

-0.4 

-0.3 

-0.2 

-0 . 1 

0.0 

+  0.1 

+0.2 

+0.3 

+0.4 

+0.5 

1  0.2361  0.737U  0.23*5  C.2394  0.2399  0.2401  0.2399  0.2394  0.2385  0.2374  0.2361 

2  0.1910  0.1978  0..'*J  +  4  0.2101  0.2  139  0.2152  0.2139  0.2101  0.2044  0.1978  0.1910 

3  0.1.39  0.1549  Q.J'.Sa  C.1758  0.1‘1/G  0.1850  0.1826  0.1758  0.1G58  0.1545  0.1439 

4  O.i:,'  'J.  1?:5  (#‘.235;  0.1403  0.1540  0.1  554  0.1540  0.1469  0.1359  0.1229  0.1109 

5  0 . 0  &  9  2  -'.1-17  0.114J  0.1249  0.130G  0.1324  0.1306  0.1249  0.1148  0.1017  0.0892 

6  0 . C  7  -« t  O.Cc/l  0.u396  O.IJCi  0.1121  0.1132  0.1121  0.1061  0.0998  0.0871  0.0744 

7  0.^:8  0.U/-.7  0.0;.-:4  0.0''49  o. 0**75  0 . 0  'j  8  2  0.C975  0.0949  0.0884  0.07G7  0.0638 

8  0..S1-;  9.-  t..»J  0.0/95  O.C+44  O.Cs'-O  0.0863  G.0860  O.OS44  0.0795  0.0GS8  0.0558 

9  C.C4‘o  C.  -Jo  0.7722  0.3753  0.0767  0.0769  0.07G7  0.0758  0.0722  0.0626  0.0496 

°  13  O.C4-o  0.057b  0.r.660  O.C--SG  C.0692  0.0693  0.0692  O.OGSG  0.0660  0.057G  0.0446 

11  0.‘J4:6  0.0534  O.CuCS  0.0070  0.0630  0.0630  0.0630  0.UG26  O.OG08  0.0534  0.0406 

12  0...372  0.M4-1)  0 . 0 > o 3  0.ft*»76  0.»*r.77  0.0578  0.0577  0.057G  0.05G3  0.0499  0.0372 

11  0.0J43  U.C401  0.0523  C.0552  0.6533  0.0533  0.0533  0.0552  0.0523  0.0469  0.0343 

14  0.031  )  9.0442  C .  C  ’« 8  3  0.C494  C.U4‘)5  0.0495  0.0495  0.0494  0.0488  0.0442  0.0319 

15  n.02“8  O.r.::  0.C457  O.O462  0.0462  0.0462  0.0462  0.04G2  0.0457  0.0418  0.0298 

16  0.0279  0.63//  0.0450  0.C453  0.0433  0.0433  0.0435  0.0433  0.0430  0.0397  0.0279 

17  0.U2C?  0.0378  0.1*4»6  0.0*468  0.0408  0.0408  0.0408  0.0408  0.0406  0.0378  0.0263 

18  0.024a  0.03*1  0.0384  l'. 0335  0.0335  0.0385  O.C3fc5  0.0385  0.0384  0.0361  0.0248 

1  »  .9255  0.0545  0.03t»4  0.0365  0.0JG5  0.0365  0.0365  0.0365  0.0364  0.0345  0.0235 

29  0.0223  0.0550  Q.054G  0.0347  0.0347  0.0347  0.0347  0.0347  0.034G  0.0330  0.0223 


Table  3  /( x)  for  d=~0.5  and  d=-0.4 

X 


0.0 

0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

1.0 

1 

0.551 

0.039 

0.734 

0.833 

0.934 

1.033 

1.123 

1.202 

1.2G2 

1.300 

1.313 

2 

0.147 

0.215 

0.312 

0.443 

0.635 

0.871 

1.1G0 

1.476 

1.775 

1.994 

2.075 

3 

0.030 

0.054 

0.098 

0.17G 

0.312 

0.545 

0.920 

1.468 

2.145 

2.759 

3.015 

4 

0.005 

0.012 

0.027 

0.059 

0.130 

0.233 

0.003 

1.221 

2.238 

3.425 

4.003 

5 

O.COO 

0.002 

0.007 

0.013 

0.049 

0.133 

0.353 

0.904 

2.100 

3.933 

5.000 

°  6 

0.000 

0.000 

0.0U2 

0.605 

0.018 

0.059 

0.194 

0.G21 

1.830 

4.270 

6.000 

7 

0.000 

G.COO 

0.000 

0.002 

o.coG 

0. 025 

0.103 

0.408 

1.513 

4.443 

7.000 

e 

C.  000 

6.000 

0.0U  J 

O.COO 

0.002 

0.011 

0.u53 

0.259 

1.204 

4.472 

8.000 

9 

0.000 

0.0  JO 

0.000 

0.000 

0.  OUO 

0.004 

0.027 

0.161 

0.932 

4.382 

9.000 

10 

0.000 

u.  000 

0.000 

O.wOQ 

0.000 

6 

0.002 

*  -0.5 

0.013 

6.099 

0.707 

4.200 

10.000 

1 

0.597 

0.G85 

0.778 

0.872 

0.9C4 

1.049 

1.122 

1.178 

1.213 

1.22G 

1.213 

2 

0.181 

0.20  3 

0.373 

0.533 

0.73b 

0.977 

1.244 

1.49G 

1.G80 

1.748 

1.610 

J 

O.OU2 

0.07G 

0.136 

0.242 

0.423 

0.714 

1.139 

1.6G5 

2.141 

2.340 

2.161 

4 

0.00 

0.019 

0.CU3 

0.094 

0.205 

0.437 

0.88  5 

1.622 

2.483 

2.902 

2.483 

5 

0.002 

0.605 

0 .  C 1 2 

0 .  C  3 '« 

0.091 

0.242 

0.618 

1.43G 

2. GOO 

3.420 

2.G90 

°  6 

0.000 

0.001 

0.004 

C.012 

0.039 

0.  12G 

0.404 

1.191 

2.778 

3.904 

2.778 

7 

C.OJO 

0.000 

O.GjO 

0.004 

O.Olu 

0.064 

0.254 

0.943 

2.7.69 

4.3G3 

2.769 

8 

0.000 

0.000 

0.039 

0.0J1 

0.636 

0.032 

0.I5G 

0.724 

2.088 

4.101 

2.686 

9 

0.003 

0.003 

0.3-0 

0.030 

0.065 

0.016 

0.094 

0.  ;43 

2.553 

5.244 

2.S53 

10 

0.000 

0.000 

0.300 

0.  Quo 

0.601 

0.068 

0.05b 

3.401 

2.384 

5.677 

2.366 

l  -  -0.4 
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Table  4  f(x)  for  3  =  -0.3  and  8=  -0.2 


X 


0.0 

0.  1 

0.2 

0.  3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

1 . 0 

1 

0.649 

0.737 

0.826 

0.913 

0.993 

1.012 

1.11b 

1.14  9 

1.161 

1.14  9 

1.116 

2 

0.231 

0.332 

0.469 

0.C45 

0.639 

1.09  3 

1.515 

1.4/7 

1.357 

1.477 

1.515 

3 

0.064 

0.119 

0.204 

0 . 59  6 

0.  002 

G .  5 1>  1 

1.404 

1 . 8  b  6 

1.973 

1.665 

1 . 4  b  6 

4 

0.016 

0.053 

0.075 

0.170 

0  ,  3  b  5 

D.  735 

1.347 

2.04.1 

2 . 4  C  9 

2.  bbl 

1.34/ 

5 

0.004 

0.910 

0.028 

0.0/6 

0.261 

0.515 

1.14  2 

2.255 

2.1  y, 

2  .  <  5  > 

.  .  .  2 

a  6 

0.000 

0.003 

0.  CIO 

0.052 

0.  1  Co 

0.339 

O.Vjo 

2.329 

3.2/  2 

2.5.  > 

6  .  J  9  o 

7 

0.000 

0.000 

0.003 

0.014 

0.033 

0.21b 

6  .  t>  o  5 

2.55/ 

5.  7  i  5 

2  .  3  7 

o .  o  v  5 

8 

0.000 

O.OGC 

0.001 

0.006 

0.628 

0.133 

6.627 

2.32  7 

*4  .  1  «.  > 

2.527 

4.8.7 

9 

0.000 

0.000 

0.000 

0.002 

G.G14 

1).  -o5 

6.4  79 

1.251 

4  .  u  5 

2.261 

C  .  4  7  j 

10 

0.000 

0.000 

0.GUQ 

0.000 

0.007 

i 

0.050 

*  -0.  3 

0. 360 

2.138 

5.092 

2.15o 

6 .  luO 

1 

0.709 

0.704 

0.373 

0.955 

1.022 

1.073 

1.105 

1.110 

1.10  5 

1.673 

1 .  C  2  2 

2 

0.304 

0.429 

0.591 

0.786 

1.001 

1.203 

1.351 

1.40b 

1.551 

1.263 

I  .  001 

3 

0.104 

0.184 

0.321 

0.542 

0.366 

1.  2G6 

1.628 

1.779 

1 .628 

1 .  2  ..  4 

0.  Cb6 

4 

0.032 

0.C71 

0.155 

0.330 

0.663 

1.22*4 

1.  S74 

2.190 

1.874 

i .  2 :  *. 

C  .  b  6  8 

a  3 

0.010 

0.020 

C.  070 

0.186 

0.4/5 

*.103 

2.066 

2.627 

2 . 6  b  6 

1.165 

0  .  *  7  5 

°  6 

0.003 

0.009 

0.030 

0.100 

0.319 

0.940 

2.194 

3.083 

2.154 

0.540 

6.519 

7 

0.000 

0.003 

C.013 

0.052 

0.267 

0.  708 

2.265 

3.553 

2.266 

0  .  7  C  o 

0.  2C7 

8 

0.000 

G.  001 

0.005 

0.027 

0.131 

0.607 

2.255 

4.6  33 

2.255 

6.6*7 

6.151 

9 

0.000 

0.000 

0.002 

0.013 

0.081 

0.468 

2.201 

4.520 

2.201 

0  .  *4  b  0 

6.031 

10 

0.000 

0.000 

0.000 

0. 0o7 

0.049 

0.354 

2.105 

5.012 

2.  165 

0.354 

6. 6*49 

6  »  -0.2 


Table  5  f(x)  for  d=  -0.1  and  d  =  0.0 


X 


0.0 

0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0,8 

0.9 

1.0 

1 

0.776 

0.858 

0.933 

0.998 

1.048 

1.080 

1.091 

1.050 

1.048 

C.99S 

C.933 

2 

0.407 

0.561 

0.747 

0.950 

1.143 

1.283 

1.335 

1.283 

1.143 

0.950 

0.  747 

3 

0.174 

0.304 

0.514 

0.820 

1.199 

1.542 

1.685 

1.542 

1.199 

0.820 

0.514 

4 

0.068 

0.148 

0.316 

0.640 

1.174 

1.796 

2.099 

..750 

1.174 

0 . 0  4  0 

0.316 

a  5 

0.025 

0.008 

0.130 

0.401 

1.072 

2.00  7 

2.552 

2.007 

1.072 

0.  «61 

0.160 

0.009 

0.030 

0.098 

0.313 

0.923 

2.154 

3.027 

2.154 

C.923 

0.313 

C .  v  5  3 

7 

0.005 

0.013 

0.052 

0.205 

0.760 

2.230 

3.514 

2.230 

0.760 

0.265 

0.062 

8 

0.001 

0.005 

0.027 

0.130 

0.603 

2.240 

4.007 

2.240 

0.603 

0.136 

0.027 

9 

0.000 

0.002 

0.013 

0.031 

0. 466 

2.193 

4.505 

2.193 

0 . 4b6 

0.081 

C.013 

10 

0.000 

0.000 

0.007 

0.049 

0.353 

6 

2.101 

«  -0.1 

5.602 

2.1C1 

0.353 

0.045 

0.007 

1 

0.851 

0.926 

0.990 

1.040 

1.071 

1.082 

1.071 

1.040 

0.990 

U.92o 

0.851 

2 

0.551 

0.734 

0.934 

1.123 

1.262 

1.313 

1.262 

1.123 

0.934 

0.734 

0.551 

3 

0.299 

0.505 

0.807 

1.179 

1.517 

1.657 

1.517 

1.173 

0.8C7 

0.505 

0.229 

4 

0.147 

0.312 

0.633 

1.160 

1.775 

2.075 

1.775 

1.160 

G.633 

0.312 

0.147 

5 

0.067 

0.179 

0.458 

1.064 

1.993 

2.534 

1.953 

1.064 

0.458 

0.179 

0.067 

°  6 

0.030 

0.098 

0.312 

0.920 

2.145 

3.015 

2.145 

0.920 

0.312 

0.09  8 

0.030 

7 

0.013 

0.051 

0.204 

0.758 

2.226 

3.506 

2.226 

0.758 

0.204 

0.651 

0.015 

0 

0. 065 

0.027 

0.130 

0.663 

7.238 

4.003 

2.238 

0.663 

0.130 

0.627 

0.605 

9 

0.002 

0.013 

0.031 

0.466 

2.192 

4.  501 

2.152 

0 . 4  66 

U.GS1 

0.015 

0.002 

10 

O.OOU 

0.007 

0.049 

0.353 

2.100 

5.000 

2.100 

0.353 

0.049 

0.007 

0.000 

*5  *  C.O 


where  /i  is  the  mean.  By  removing  the  absolute  sign  from  Equation  (8),  r  can  be 
as  the  sum  of  the  two  integrals, 

Hot,  =  ^  (/«-a :)f(x)dx  +  \  (.v-/i)/(-v)</.v. 

Nothing  that  is  independent  of  the  variable  .v  and  using  relationships  such  as 

\’/U)rfjr-l-F(/i) 


and 


written 

(9) 
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\  X  /(,)</,  ft  -  \  \f(x)l/x. 

Equation  (9)  becomes 

rUt.  o)  -2’if  (;t)-2\  ‘  x/(x)dx.  (10) 

Denote  the  integral  in  Equation  (10)  by  //(,««)  and  for  ft">0, 

,/(,)</, 

=  L'vC(ft’  5)cosh-M-v+5-0.5)]  dX 

=  G(rt,  «)\o  coshi  [a(.v  +  5— 0.5)]]  dx' 

The  integral  is  similar  to  the  one  in  the  derivation  of  the  mean  with  one  exception,  the 
upper  limit  is  now  /t  instead  to  1.0. 

Evaluating  this  integral  with  proper  limits,  H(/j)  becomes 

*  - “"h  c*«-o.5+*n+io.  £„“schi(s-o.5+°n 

for  «>0  and  becomes 

//(,,)=  p,</,=0.125  =  0.125  for  «  =  0. 


Therefore  H(;i)  is 


Eig.  I  Equi  mean  curves  (dotted  lines),  and 
cqui-absolute  deviation  curves  (solid 
lines).  The  numbers  in  parenthesis 
are  the  values  of  normalized  absolute 
deviation,  y 


=  G(a\  S}  a,i  tanh  [>(5-0.5  +  //)] 


H(,t) 


4-lr.a  c°sh  1X5— 0.5)] 
g  cosh  [X5-0.5  +  /4] 
for  or>0 


=  0.125  for  a  —  0. 


The  absolute  deviation,  r.  for  the 
FRPDF  is  therefore 


r(«,  8)  =  2,tF(n)-2H(n)  (12) 

where  the  explicit  expressions  for  ft,  F{ft), 
and  H{fi)  are  given  in  Equations  (7),  (4) 
and  (11)  respectively. 

By  use  of  Equations  (7),  (12)  and  the 
aid  of  the  computer,  the  mean  and  the 
absolute  deviation  are  computed  for  a  =  0, 
1,  •••,  20  arttt-^=-0.5,  -0.4,  •••,  +0.4, 
+  0.5  and  the  equi-mean  and  equi-absolute 
deviation  curves  for  a  =  0,  1,  •••,  10  are 
given  in  the  nomograph  in  Fig.  I. 

By  computing  the  mean  and  the  ab¬ 
solute  deviation  of  sample  data,  one  can 
find  the  values  of  a  and  5  directly  from 
the  nomograph.  Upon  substituting  these 
values  of  at  and  8  into  Equation  (4)  one 
will  obtain  "the  best  fitted"  probability 
distribution  function  for  this  particular 
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set  of  samples  according  to  the  criterion  cited  earlier  (/>.,  this  set  of  samples  has  the  same 
mean  and  the  same  absolute  deviation  as  that  of  the  sample  having  the  probability  distribu¬ 
tion  function  with  the  values  of  a  and  3  taken  from  the  nomograph). 


6.  The  Variance  and  Moment  Generating  Function 

From  a  practical  viewpoint  the  variance  and  moment  generating  function  are  of  little 
use;  however,  they  arc  investigated  for  theoretical  considerations  to  verify  that  the  variance 
is  finite  and  that  the  moment  generating  function  exists. 

By  definition  the  variance,  a2,  is: 

<F=\  ( x—/t)2f(x)dx .  (13) 

Equation  (13)  ca..  be  simplified  as  follows; 


2=  \  x2/(x)dx :-/r. 


(14) 


When  <ar  =  0,  /(*)=  1,  0<x£;l  and  /i=0.5,  therefore, 

,  1 
12  * 

This  is  correct  for  the  uniform  distribution. 

For  or>0,  the  integral  in  Equation  (14)  is  difficult  to  evaluate  and  since  it  is  shown 
that  r  is  used  in  lieu  of  a2,  with  the  nomograph,  for  finding  parameters  a  and  d,  only  the 
existence  of  the  variance  will  be  shown.  Then  from  Equation  (14) 


\'  x2f(x)dx-,r  =  C(a,  S)\0 


r*  r  i  f*'o+o.5i 

r 

J a 


cosh2  [>(*  +  5-0.5)] 
dz+  2(0  5~S) 

ft-  .'o 
]■ 


dx  —  (C-. 


cosh2  c 


dz 


cosh2 c 

+  (0.5-5)2  1 

Oi  L'i-o.r.1  coshJz  1 

The  first  integrant  can  easily  be  shown  to  be  integrable  by  use  of  the  theorem:  A  func¬ 
tion  continuous  on  a  closed  interval  is  integrable  there-6!!.  Also,  since  it  is  a  proper  distribu¬ 
tion  function,  it  has  a  unique  characteristic  function,  by  Levy's  Theorem. 

Let  us  denote  the  moment  generating  function  for  .v  by  MAO).  By  definition, 

MA0)=\  f{x)enidx  (15) 

When  n  =  0,  /(*)=!,  0<*<1,  therefore 

MA0)  =  \  'e"'dx^  e"~]  . 

When  a>0.  Equation  (15)  becomes, 

MA0)^\'c(a,  3)  .  ,  -e  'd\ 

;l>  cosh-poLr  +  o  —  O.o)] 

G(«,  cosh- 1  r>(.v+  o  0.5)1  </A' 

Using  the  same  argument  of  continuity  and  integrability  of  this  integrant,  it  follows  that  the 
moment  generating  function  MAO)  exists. 

To  facilitate  application  of  the  FRPDF,  the  following  tables  and  a  .Nomogram  is  pro¬ 
vided.  Also,  several  figures  showing  the  shape  of  the  FRPDF  w  ith  different  values  of  param¬ 
eters  a  and  3  are  given.  The  attached  Appendix  will  help  the  reader  get  a  first-hand  fee* 
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of  its  ease  in  application  and  goodness-of-fit. 


7.  Conclusions 

The  objective  of  this  paper  has  been  to  present  a  finite  range  probability  distribution 
function  which  was  developed  by  Braswell  and  further  refined  and  applied  by  Manders. 
The  probability  density  function,  /(.v),  is  strictly  uniniodal  as  shown  in  Figs.  2  through  7. 
When  S  0,  the  probability  density  function  is  always  symmetrical  with  respect  to  the  vertical 
line  .v  —  0.5,  and  the  distribution  function  is  a  family  of  S-shaped  curves.  When  5<0,  the 
peak  of  the  probability  density  function  shifts  toward  the  x  =  1  line  and  when  5>0,  it  shifts 
to  the  .v  —  0  line.  When  a-=0,  the  FRPDF  distribution  is  identically  uniform;  for  the  large 
value  of  a,  a=10  or  more,  the  peak  of  the  distribution  becomes  narrower  and  higher  at 
appropriate  value  of  .«  for  given  value  of  3.  Eventually,  this  peak  will  reach  infinity  in  the 
limit,  and  the  function  becomes  an  impulse  function  at  a  certain  value  of  x.  This  is  very  con¬ 
venient  in  practical  cases.  Since  this  indicates  that  for  a  certain  statistical  sample,  if  the  value 
of  n  is  large,  this  statistical  sample  can  be  treated  approximately  as  a  deterministic  one. 
Flencc  3  can  be  considered  as  the  location  parameter,  while  a  can  be  considered  as  the 
shape  parameter. 

The  new  FRPDF  can  be  easily  applied  to  almost  any  practical  problems  where  experi¬ 
mental  data  are  easily  collected.  Once  the  necessary  data  are  tabulated,  one  can  easily  compute 
the  mean  and  the  absolute  deviation  of  that  sample  data.  Then  by  use  of  the  Nomogram  in 
Fig.  1  the  values  of  the  parameters,  a  and  S,  can  be  found  directly.  Substituting  these 


Fig.  2  fix)  for  the  FRPDF  <J-0.0  Fig.  3  f(x )  for  the  FRPDF  «J=-0.1 

If  d=  fO.J,  replace  x  by  l-x 
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values  of  it  and  d  into  liquation  (4)  one  can  obtain  “the  best  fitted"  probability  distribution 
function. 
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ON  TESTING  AND  APPLICATION  OF  A  NEW  FINITE  RANGE 
PROBABILITY  DISTRIBUTION  FUNCTION 


By  Robert  N.  BkAswti.i. 

University  of  Florida 
Gainesville,  Florida 

Clara  Fu-Mei  Mandlrs 
Radiation,  Incorporated  Melbourne,  Florida 


To  demonstrate  flexibility  and  sensitivity  this  paper  will  cover  some  aspects  of  testing 
and  application.  There  have  been  several  other  experiments  using  the  new1  FRPDF  and  in 
all  cases  it  has  proven  superior  in  results  and  much  easier  to  use.  Even  when  the  distribu¬ 
tion  of  the  data  was  known  the  FRPDF  was  preferable.  The  complex  expression  did  not 
deter  new  users  from  being  attracted  to  further  applications. 


1.  Testing  and  Application 


An  important  problem  in  statistics  is  to  find  how  well  a  sample  taken  from  a  population 
agrees  with  some  distribution  function  assumed  for  that  population.  Two  such  tests  are 
considered  here. 


The  range  of  .v  is  divided  into  M  equal  regions  and  the  number  of  sample  points  falling 
within  each  region  is  counted.  Let  Yu  Yz,  •••,  be  the  result.  From  the  assumed  distribu¬ 
tion  and  the  size  of  the  sample,  the  expected  number  of  points  in  each  region  is  computed: 
gi,  gt,  •,  gt.t.  The  deviation  between  this  and  the  actual  result  is  expressed  by 

d=Y  (D 

Li  £»» 


■where 


^H1  fa  P  A I 

=  the  number  of  expected  points  in  the  w-th  interval 
.w 

m  1 

=  sample  size 

y„,  =  the  number  of  sample  points  in  the  w-th  interval 
Pm  =  \  "  f(x)dx  =  F(x„)  -  F(.r„,  , )  /«  =  1 ,  •  •  • ,  M 

-  xm  ] 

Ffam,)-Q  for  iu=  1 

the  probability  of  sample  points  falling  in  the  m-th  interval. 
This  deviation  is  used  to  ascertain  the  confidence  level  of  the  assumed  distribution. 


*  Rcccocd  24.  l-cb.,  IT70. 

(R apfuintzd  fa om  Rep.  Stat.  A ppl.  Rte.,  JUSE,  l lot  17,  No.  2,  1970,  ba 
peluruAbion  of,  the  pubLL&heA. ) 
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As  an  application  of  the  method  described  in  this  paper,  consider  a  class  grade  distribu¬ 
tion  of  50  graduate  engineering  students.  Table  1  gives  the  grades  for  the  students  [1], 

Srip  1  Compute  an  estimate  of  the  normalized  mean,  /i,  and  the  absolute  deviation, 
from  the  data  given  in  Table  1.  For  convenience  in  analysis  we  will  normalize  the  grade 
range  from  [0,  I00J  to  CO,  1].  They  are  found  to  be; 

/« =  0.829 
r  =0.093. 


Table  1  Class  grades  of  a  graduate  engineering  class  of  SO  students 


< 

Gr 

/ 

Gr. 

i" 

1 

81 

11 

85 

21 

2 

86 

12 

69 

22 

3 

78 

13 

95 

23 

4 

85 

14 

87 

24 

5 

79 

15 

69 

25 

6 

85 

16 

85 

26 

7 

92 

17 

90 

27 

8 

83 

18 

71 

28 

9 

96 

19 

68 

29 

10 

89 

20 

79 

30 

Step 

II  Locate 

the 

intersection  of 

estimate  the  coordinates  of  this  intersection. 


Gr. 

i 

Gr. 

1 

i 

Gr. 

80 

31 

67 

41 

78 

49 

32 

89 

42 

95 

96 

33 

77 

43 

86 

63 

34 

95 

44 

75 

92 

35 

63 

45 

96 

66 

36 

100 

46 

96 

52 

37 

85 

47 

98 

79 

38 

97 

48 

83 

94 

39 

85 

49 

81 

96 

40 

90 

50 

91 

i 


829  curve  with  f  =  0.093  curve  on  Fig.  1,  and 


Fig.  I  The  FRPDF  F(jr)  with  a  =  5.5, 
i~  -0.4.  Fitted  to  curve  of 
class  grade 


This  is  equivalent  to  the  estimation  of  two 
parameters  a  and  8  in  our  F(x).  Here 
«-=5.5 
8  *  —0.4. 


Step  III  Substituting  these  values  into 
Equation  (4),  we  obtain 

.  cosh (0.55) sinh (5.5*)  ... 

W  sinh  (5.5)cosh  C5.5(jr— 0.9)]  '  ' 1  ’ 


Step  IV  Divide  the  range  of  x  into  10 
equal  reg!  >ns:  M  =  10 

1)  Compute  F(xm): 

\  =  cosl1  sinh  (5.5*.,) 

KXm)  sinh  (5.5)  cosh  5.5(*- -0.9) 

m=  1,  •••,  10. 

With  the  aid  of  the  computer,  the  values 
of  F(xJ  are,  0.0001,  0.0005,  0.0017,  0.0054, 
0.0161,  0.0473,  0.1329,  0.3328,  0.6664,  and 
1.000. 

2)  Compute  pm\ 

Pm  =  \  '  f(x)dx  =  F(xm)~  F(xm-x), 


i 
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m  I,  2,  M. 

Using  the  values  for  /  (a.,),  the  values  ol  l‘  are,  0.0001,  0.0004,  0.0012,  0.0037,  0.0107, 
0.0312,  0.0856,  0.1999,  0.3336,  and  0.3336. 


Siii-  V  Prepare  Table  2. 

t  able  2  Distribution  of  class  grades 


m 

1  m 

v  Sp„ 

Km 

1 

0 

0.005 

0.005 

-» 

0 

0.020 

0.020 

3 

0 

0.060 

0.060 

4 

0 

0. 185 

0. 185 

S 

1 

0.535 

0.404 

6 

1 

1.560 

0.201 

m 

.I. 

V’.,  \pn 

(  I'm  l  , 

I'm 

7 

7 

4.280 

1 . 729 

X 

9 

9.  995 

0.099 

9 

17 

16.680 

0.0J6 

10 

15 

16.680 

0.  i69 

Total 

50 

50 

2.878 

2.  The  Z^-Test  [2,  3] 


The  purpose  of  the  test  is  to  find  the  probability  that  the  observed  deviation  between 
the  theoretical  model  and  the  empirical  data  is  in  fact  due  to  the  random  nature  of  this  set 
of  data. 

If  the  y„  are  sufficiently  large,  say  more  than  10,  the  deviation  we  defined  in  (1),  is 
distributed  according  to  the  Z-  distribution  w  ith  M-  1  degrees  of  freedom.  This  is  K.  Pear¬ 
son's  test  function  which  gives  great  weight  to  those  deviation  squares  (y„  —  #„)-  that  cor¬ 
respond  to  small  /?,.  If  the  assumed  parent  distribution  is  not  completely  known  and  k 
parameters  defining  it  have  been  determined  to  fit  the  sample,  the  number  of  degrees  of 
freedom  is  reduced  to  M—  1—  k. 

Now  we  will  apply  the  Z  -test  to  our  example. 


u 


=  2.878 

and  the  degrees  of  freedom  are  7  (10-1-2). 

For  7  degrees  of  freedom,  this  deviation  is  exceeded  about  90  percent  of  the  time  '. 
The  assumption  of  our  FRPDF  is  therefore  very  good.  There  is  thus  nothing  in  the  value 
of  X-  to  lead  us  to  reject  our  hypothesis. 


3.  The  Kolmogorov-Smirnov  Test  f3J 

It  is  also  desired  to  investigate  how  well  our  empirical  data  tits  our  theoretical  distribu¬ 
tion  by  A  -  .v  test.  The  k  —  \  test  allows  us  to  place  confidence  level  on  the  positive  as  well 
as  negative  deviations,  /.<■.,  it  allows  us  to  check  the  theoretical  distribution  for  points  of 
excessive  as  well  as  inadequate  probability. 

Define  the  one-sided  deviations  as: 

£).v*  =  sup  CF(.v„)-5„(.v„)j  and  D  sup  [£.(*.,)- F( a.)]  (4) 

Oil  v  I  1 

where 
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S...(v,„)  :  ^  and  F(.v„) 

arc  given  in  Step  IV. 

According  to  Smirnov's  asymptotic  distribution3, 

PI)  '£}'*•  '■  and  P[DS  ,'X)'e  (5) 

Now  \sc  proceed  to  test  the  example. 

D s  •  .  .0  0.0014  and  Ds  a  0.0471. 

lable  3  Calculation  of  the  one-sided  deviations,  £).>'  and  D.\~ 


m 

/  ( v„) 

i'ut.V-) 

l)s 

i). 

m 

Firm) 

Sm(Vm) 

Os' 

/>v- 

1 

0.0001 

0.00 

6 

0.0473 

0.04 

0.0073 

7 

0  0005 

0.00 

7 

0. 1320 

0. 18 

0.0471 

3 

0.0017 

0.00 

8 

0.  3328 

0.36 

4 

0.0054 

0.00 

9 

0.6664 

0.70 

5 

0.0161 

0.02 

10 

1.0000 

1.00 

See  lag.  I  tor  comparisons. 


4.  The  Glass  Bottom  Boat  Problem 

Consider  the  "glass  bottom  boat  problem"  which  is  similar  to  the  well  known  “Newsboy 
Problem."  1 

Snili’incm  of  the  Problem : 

A  photographer  must  decide  how  many  pictures  to  print  each  time  after  he  takes  a 
picture  of  the  passengers  in  the  glass  bottom  boat  at  Silver  Springs,  Florida.  Suppose  that 
the  total  number  of  the  passengers  is  fixed,  say  forty.  The  <  ost  of  a  picture  is  C  and  the 
selling  price  is  S.  An>  pictures  not  sold  at  the  end  of  the  day  are  a  total  loss.  Let  ply) 
be  the  probability  that  i  pictures  will  be  demanded  each  time.  Then  his  expected  profit  for 
each  boat  if  he  prints  li  pictures  is 

t:\  Plh)}  sSlyi*)  )  r  Shi  V  p(  i  )J  —  C/i 

v  V  n  •  1 

since  the  revenue  received  is  Sy  if.i  h,  and  is  Sh  if  y>h.  The  problem  is  to  dele-mine 
the  value  of  h  which  maximizes  his  expected  profit. 

Often  it  is  convenient  to  treat  h  and  the  demand  variable  y  as  continuous.  Then  if 
fly)  is  the  density  function  for  demand  and  Fly)  is  its  distribution  function,  the  expected 
profit  for  each  boat  when  h  units  are  printed  is 

/. [  /’( /;) J  S\  t Jl\)dy  +  Sh\  f!y)dy-Ch.  (6) 

The  optimal  h  is  then  a  solution  to  i/F^PIh)]  dh= 0. 

Lsing  l  eibnit/'s  Rule,  we  obtain 

dl':  -  o =s-c~ suit). 

dn 

Thus  the  optimal  h%  satisfies  the  equation 

FUi)-S~C.  (7) 
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Ld  S  SI .20  and  (  $0.50.  Then  /•<//)  7  12.  or  0.58.1. 

I  quation  (6)  is  a  strictly  concave  function  of  A.  This  implies  that  any  relative  max¬ 
imum  of  /  |  P  |  is  the  absolute  maximum  ttnd  the  absolute  maximum  is  unique. 

Suppose  the  demands  are  normally  distributed  then  I  qualion  (7)  becomes 

•iV)  v 

where 


i  is  the  sample  mean,  and 
(7.  is  the  standard  deviation. 


If  the  demands  follow  the  bRPDI-  then  liquation  (B.7)  becomes 

....  cosh  |  o  (or  0.5)  |  sinh(zi/;)  .S'  ( 

1  sinh(a)  cosh  |  f‘(/t  •  it  0.5)1  .$ 

where  h,  is  normalized  h. 

In  order  to  determine  the  demand  the  photographer  performs  the  following  experiment. 
Tor  each  of  ten  successive  boat  loads  he  prints  forty  pictures  (the  maximum  possible  demand). 
He  then  records  the  number  of  prints  that  he  sells  to  each  boat  load  of  passengers.  In  this 
manner  he  obtains  the  following  table: 


lahle  4  Number  of  pictures  sold  to  each  boatload  of  people 


Boat  Number 
/ 

Number  of  Demands 

is 

Boat  Number 
i 

Number  of  Demands 

Vj 

1 

3fi 

6 

30 

s 

34 

7 

25 

3 

39 

8 

37 

4 

20 

9 

23 

5 

17 

10 

33 

Table  5  Calculation  of  i,  y and  ,'.17-  .1  )-' 


• 

1  : 

t  a  -  i 

(y. -  .i): 

/ 

i. 

3  ■ 

(r>  .»•)- 

1 

36 

6. 6 

43.  56 

7 

25 

4.4 

19.36 

s 

34 

4. 6 

21. 16 

8 

37 

7.6 

57.  76 

3 

39 

9.6 

92.  16 

9 

23 

6.4 

40.96 

4 

20 

9.4 

88.  36 

to 

33 

.3.6 

12.96 

5 

17 

12.4 

153.76 

6 

30 

0.6 

.36 

Total 

294 

65.2 

530. 40 

We  find 

that  the 

estimator 

for  /(,  is  the 

sample  mean 

i , 

The  estimator  for  n,2  is: 


and  that  for  a.  is. 
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The  estimator  for  y,  is  the  sample  absolute  deviation  \y—y\, 

10 

l.~|  ;oYj  V.->|  =  6.52. 

I  1 

Now  we  are  in  the  position  to  calculate  h'  by  assuming  that  F(/t)  in  Equation  (7)  is  1) 
normal  distribution  function  and  2)  FRPDF. 

1)  Assume  F(h)  is  normal  distribution  function.  Then  h*  is  the  solution  to  Equation 
(8),  or 

•c )-«•>■ 

Hence  from  the  normal  tables 

A  —  29.4  ...  . 

7.68  =0'21  °r  /,=  3L 

Thus  31  pictures  should  be  printed. 

2)  Assume  F(h)  is  FRPDF.  In  this  model  the  variable  is  normalized;  we  should  make 
a  litisr  transformation  on  each  sample  value. 

We  find  the  normalize  y,  ix  to  be: 

-to  4 

and  the  normalized  f  to^be 


From  Fig.  I,  we  can  estimate  the  parameters  of  our  FRPDF  a  and  d.  They  are  found 
to  be 

d  =  3 
§=—  0.4. 


Table  6  Calculation  of  the  expected  profit  for  three  different  decisions 


Number  of  Pictures  Sold 

>'» 

*=>■  =  29 

h  =  31  (Normal) 

|  h  =  33  (FDPDF) 

1 

36 

29 

31 

33 

2 

34 

29 

31 

33 

3 

39 

29 

31 

33 

4 

20 

20 

20 

20 

5 

17 

17 

I’ 

17 

6 

30 

29 

30 

30 

7 

25 

25 

25 

25 

8 

37 

29 

31 

33 

9 

23 

23 

23 

23 

10 

33 

29 

3* 

33 

Total  sold 

259 

270 

280 

Total  printed 

290 

310 

330 

Total  return 

S3 10. 80 

S324.00 

S336.00 

Total  cost 

145.00 

155.00 

165.00 

Total  profit 

165.80 

169.00 

171.00 
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Substituting  these  values  into  Equation  (9)  and  rewriting  it,  we  have 
coth(jA,)—  — tunh(  — 2.7)+  ',2  cos„”!?™  "h(3, 
=  0.99101+  12  ( l  .0453) 


7  (7.4735X10.018) 


=  1.01494 


or 


tanh(3/iu)  =  0.98528 
3A0  =  2.46 


2.46 

3 


and 


->  46 

A=  3  ✓.  40  33. 


Thus  33  pictures  should  be  printed  in  this  case. 

Of  the  three  methods  employed  we  see  that  the  new  FRPDF  gives  the  BEST  decision 
as  to  number  of  pictures  to  be  printed.  More  involved  experiments  with  the  Newsboy 
Problem,  the  Glass  Bottom  Boat  Problem,  etc.,  yield  comparable  results  f6j. 
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(The  reverse  of  this  •'  i  ■-  ;s  blank) 


